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Abstract 



With the eigen-functional bosonization method, we study one-dimensional 
strongly correlated electron systems with large momentum (Ikp and/or 4:kp) 
transfer term(s), and demonstrate that this kind of problems ends in to solve 
the Eikonal-type equations, and these equations are universal, and indepen- 
dent of whether or not the system is integrable. In contrast to usual pertur- 
bation theory, this method is valid not only for weak electron interaction, but 
also for strong electron interaction. Comparing with exact solution of some 
integrable models, it can give correct results in one-loop approximation. This 
method can also be used to study electron-phonon interaction systems, and 
two coupled spin chain or quantum wire systems. 
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I. INTRODUCTION 



In contrast to two- and three-dimensional electron systems, one-dimensional electron 
system has a prominent feature that its Fermi "surface" is composed of two points ±kp, 
defined by the Fermi momentum kp, and its Hilbert space is strongly suppressed. For a 
one-dimensional electron gas, there only exist two kinds of electron excitations: one is near 
these two Fermi levels ^kf, respectively, usually called as small momentum excitation, and 
another one is the excitation between these two Fermi levels, called as large (crossover) 
momentum excitation. If there is only the former excitation, the system generally shows 
the Luttinger liquid behavior. After turning on the latter excitation, the system becomes 
an insulator for repulsive electron interaction. 

In low energy regime, the physics property of a one-dimensional electron system is de- 
termined by the electron states near its two Fermi levels ^kp. Therefore, it is reasonable 
to linearize electron energy spectrum near these two Fermi levels IQ-Q with a band-width 
parameter D, and the electrons are divided into two parts: one presents the electrons near 
the Fermi level +A;_f, called right-moving electrons with the energy spectrum eR{k) = vpk, 
and another one presents the electrons near the Fermi level —kp, called left-moving electrons 
with the energy spectrum eilk) = —vpk, where vp is the electron Fermi velocity. It is shown 
that the high order terms like /c", n =2,3,..., only contribute high order irrelevant terms 
to the effective action of the system, and can be neglected. The crossover (large momen- 
tum) excitation corresponds to the excitation between the rigth- and left-moving electrons 
with 2kp and/or 4kp (or higher) momentum transfer. In some special cases, the crossover 
excitation is very important, and determines the low energy behavior of the system, because 
the crossover excitation term is relevant, and strongly affected by the electron interaction. 
Generally, usual perturbation methods fail to treating the crossover excitation term for the 
strong electron correlation. 

The bosonization method is very powerful in studying of one-dimensional interacting elec- 
tron systems, such as quantum wires, spin chains and impurity scattering in one-dimensional 
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fermion systems, where it can exactly treat electron's density-density interactions , such 
as the interaction V {x — y) pji^{x) pi^iy) . However, for example, if there meantime appear an 
interaction term V pBa{x)pL(T{x) and a crossover term A['?/'|j^(x)^/'Lo-(x) + ^/']^^(x)'?/'ij,o-(a;)] in the 
Hamiltonian of the system, it is hard to use the usual bosonization and functional bosoniza- 
tion methods to effectively treat these two terms as a whole, because the A-term mixes the 
right- and left-moving electrons which induces a highly relevant term, in the bosonization 
representation of the electron fields, it reduces to some cosine- and/or sine- terms. The usual 
bosonization methods can only tell us which ones are relevant or irrelevant, but cannot be 
used to effectively treat these terms, such as calculating spin and charge collective excitation 
spectrums, correlation functions, effective action of the system, and so on. Therefore, the 
usual bosonization methods do not provide more useful informations than usual perturba- 
tion methods. It is a big challenge to effectively treat these crossover excitation terms which 
often appear in one-dimensional strongly correlated electron systems, and desires to find a 
new method which can be used to exactly and effectively treat them. 

In this paper, we introduce a new method 0, called the eigen-functional bosonization 
method, which is very simple and effective to treating these crossover terms, and show that 
the crossover terms can be exactly treated by solving the Eikonal-type equations. It is 
wonderful that the problems of one-dimensional strongly correlated electron systems end in 
to solve the Eikonal-type equations that can be exactly treated by a series expansion and/or 
computer calculations. This method is universal, and independent of whether or not the 
system is integrable. For example, it can also be used to treat electron-phonon interaction 
systems, and two- and three-dimensional electron systems 0. In Sect. II, for simplicity, we 
introduce the eigen-functional bosonization method by studying a one-dimensional spinless 
fermion system with a staggered chemical potential or a dimerized lattice sites, and show 
that this kind of problem ends in to solve simple Eikonal-type equations. At one-loop 
approximation, our calculation is consistent with the exact solution of the quantum sine- 
Gordon model. In Sect. Ill, with this method, we study the Hubbard model at half filling, 
and calculate charge collective excitation gap which is agreement with the exact solution of 
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the Hubbard model at half filling. In Sect. IV, we study the Hubbard model at half filling 
with a staggered chemical potential, and show that there is a phase transition from band 
insulator to Mott-type insulator at the electron interaction strength Vt where the charge 
collective excitation gap is zero, and a spontaneous lattice dimerization takes place. This 
property of the system can be used to explain the origin of the ferroelectricity of transition 
metal oxides P]. The conclusion and discussions are given in Sect. V. 



II. EIGEN-FUNCTIONAL BOSONIZATION AND ITS APPLICATION 

In this section, we introduce the eigen-functional bosonization method by studying a 
one-dimensional spinless fermion system with a staggered chemical potential or a dimerized 
lattice sites. It is well-known that in usual bosonization representation of fermion fields, 
the system reduces into the quantum sine-Gordon model (QSGM) which is integrable, and 
has been exactly solved [Q. The exact solution of the QSGM is a criterion for this new 
method. We show that in one-loop approximation this method can give correct results for 
any fermion interaction strength, whereas usual perturbation methods are invalid for strong 
fermion interaction. 

Generally, we consider a strongly correlated spinless fermion system with the Hamiltonian 



H = -ifivp J rfa:[V'lj(x)a^.V'i?(a:) - V'i(a;)5a'V'L(a^)] 

dx[\ij\i{x)^L{x) + \*^i{x)tlJn{x)] +V f dxpR{x)pL{x) (1) 



where iPr(l){x) is the right (left)-moving fermion field, and Pr{l){x) = '^R(L)i^)'^R{L){x) is 
the right (left)-moving fermion density. The A-term represents a staggered chemical potential 
(A = A* = Ao), or a lattice dimerization parameter (A = —A* = iu). Generally, we can write 
A as A = Ao+m. We can also study the influence of quantum fluctuation of electron-acoustic- 
phonon interaction on the lattice dimerization by adding the term, Ua{x)[pR{x) + Pl{x)]. 
In bosonization representation, 4>R(L)ix) = {^^y^'^exp{—i^R(L)ix)}, where dx^R{L)ix) = 
±27rp^(2,)(x), a ~ 1/-D, and D is the band-width, we have the standard QSGM |TD| {u = 0), 



4 



HsG = ^j dx{[dMx)f + ^\^) + ^ cos(/3$(a;))} 

where, mV/?^ = A/(7ra/it;,), /5V(47r) = (1 - lf'^/{l + lY'^ = g, v, = vpVT^, 
7 = V/{2rrnvF), = [^r{x) - n(x) = (3[d,^R{x) + 9,<l>i(x)]/(47r), and 

[n(a;), $(y)] = ^^(x — y). The parameter m is related to the mass of the sohtons and anti- 
sohtons or breathers |]10|]. According to the relation (3"^ = Airg, at = 47r, the Hamiltonian 
(|I|) describes a free massive fermion system {V = 0). As > Att {g > 1), there is an 
attractive interaction between the fermions '4)r{x) and i/jlIx). However, due to the crossover 
term (A-term) opens a gap in fermion excitation spectrum, the system is stable. As jS"^ < Air 
{g < 1), the fermion interaction is repulsive. In the fermionization representation it 
is very clear that in the regime < (3"^ < An the QSGM represents different physics from 
that in the regime Air < (3'^ < Sn. The critical point = 47r is a phase transition point 
which separates one phase from another one. 

After introducing auxiliary boson fields (j)R{x,t) and (j)L{x,t) which decouple the four- 
fermion interaction by adding the constraints, PR{L){x,t) = il)^^^j^-^{x,t)il)R(^^{x,t), we only 
have the quadratic terms of the fermion fields ipRiL){x,t)- Integrating out them, we obtain 
an effective action, 



S[p, 0] = W[(f)] - / dtdx{(f)R{x, t)pR{x, t) + t)pL{x, t) + Vpr{x, t)pL{x, t)} 



W[(j)] = iTrln{Mo) + i / d^ / dtdx tr U{x,t)G{x,t; x' ,t' ,[^(f)]) 



I t'^t 



(2) 



where M = Mq + 0, Vr^q = idt ± ihvpd^, MG{x, t; x', t' 

I 



A*, Vl 



(j)R{x,t), 
0, 4>L{x,t) 



6{x — x')6(t — t'), and 
\ 



In order to calculate the Green's functional G{x,t; x' ,t' , [0]), we solve the eigen-functional 



equation |]13[ of the operator M, 



M^ti^,t, 



(3) 



With the orthogonality and normalization of the eigen-functionals. 
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the Green's functional can be written as, 



G{x, t; x' , t' 



EE 



(4) 



i k,uj 

If we can exactly solve the eigen-functional equation @, generally we can calculate the 
Green's functional (^), and the effective action which only includes the boson fields. 
This process is called the eigen-functional bosonization. 

At (f){x,t) = (interacting- free case), the system becomes a band insulator, and the 
eigen-equation (|]) can be exactly solved. It has two sets of solutions corresponding to the 



eigen-values -^^^''[0] = a; =F E^, respectively, where E^ = J k'^ + |Ap (choosing h = vp = I) 



<Ha:,t,[0]) = (^)V2 



-VkC 



-ie 



1 



V 



ikx—iujt 



I 



-ie 



Akx—iujt 



(5) 



where ui = 1 



'1 + k/Ek)/2, tan(2^) = u/Aq, and T and L are the time and space 



lengths of the system, respectively. 

At A = 0, the Hamiltonian (|1|) represents the Luttinger liquid, and the eigen-functional 
equation (|^) has the exact solutions with the eigen-values Ej^^ = =f \k\, respectively. 



'^ki\^,t, [0]) 



1 



QO {x,t) 



Akx—iujt 



(6) 



where Pij(2,)Q5j(L)(^? ^) — ~4>R{L){x,t), ant 9{x) is a step function, 9{x) = 1 for a; > 0, and 
9{x) = for x < 0. With the exact solution (|), we can easily obtain the effective action of 
the system, and the fermion Green's function which are the same as that in Ref. 0,0. 
In general, for 0(x, t) ^ 0, we can choose the following formal solutions with the eigen- 



values E, 



^^■''^0] = 00 ^ Ek + Ti^i^\(f)], respectively. 



TL' 



( 



TL' 



,jA;x-i(Lj+E^+'[0])t 



(7) 



where s[,'^-'[0] = Jq J dtdx'^^^^\x,t, t)\l/^'^(x, t, [^0]) are independent of 00, and 

\Ak\ ~ 1 is a normalization constant. The boson fields (5_R(L)(a;, t) and Qii(L){x,t) satisfy the 
differential equations, 



{VL + (t)L{x,t))e 



QLix,t) 



0, 



Vk 



(8) 



With the following definitions 



the differential equation (H) can be re-written as, 



{VlVr + lEkdt + /i?^ + VrJrVlIr = - ^ + D^/h PlQ° 



+ lEkdt + + VrUVlU = 

{VRVL-tEkdt + tkd,)fR + VRfRVjR= 
{VRVL-iEkdt + zkd,)h + VRfLVLh= - 



Uk 

\Muk 

Vk 

\Muk 

Vk 

\Muk 

Uk 



(9) 



where Q^{x,t) = Q^R{x,t) — Ql{x,t). This is usual Eikonal-type equations [|1^,|T^], and can 
be solved by a series expansion of Q^{x,t), or by computer calculations. These equations 
are universal, and independent of whether or not the system is integrable. For example, 
after including the electron-acoustic-phonon interaction, they keep invariance if one replaces 

0R(L)(x,t) by 0ij(L)(x,t) - Uaix,t). 



The non-linear terms will produce the high-order terms of {\\\/hvF)Q^{x^t)^ and con- 
tribute cubic or higher order terms of the boson fields t) to the effective action, as 
\X\/hvF ^ 1 they can be neglected. Only keeping the linear terms (one-loop approximation), 
we can obtain 



fL{q,n) 



V 



\ 



1 



( 



n^-q^ + 2{EkVt - kq) 



^2_q2_ 2{EkVL + kq) 



^iSi + q) 



Q{q,n) 



\M'"k 



(n-q) 



Q{q,n) 



(10) 



where Q{q, f2) = 4>R{q, fi)/(fi — g) — 4>L{q, ^)/{^ + q). Generally, these boson fields depend 
on the fermion momentum k. However, the low energy physical property is determined by 
the fermions near the Fermi levels ±kp, we therefore can take the approximation in (|10|), 
~ and Ek ~ |A|. Under this approximation, these boson fields are independent of the 
fermion momentum k. 

According to equations (Q) and (|^), the Green's functional can be written as. 



G{x,t;x',t, [0]) 



Gr(x - x')e««(^'*)-««(^''*), Zix - x')e««(^'*)-«^(^''*) 



V 



:ii) 



Z{x - x')e^L{x,t)-Qn{x',t)^ ^^^^ _ ^yQU-,t)-QUx',t) 

where GniL){x) = G^^^^ix) ± gix), G^^^^ix) = ±l/(27r2x), gix) = lSlndk{^^ - 

and Z{x) = ^ rffc-^e*^^, where D is the band-width. Here we have used QR(L)(a;,t) = 

~QR{L){x,t), at one-loop approximation this relation is exact. With the relation 

1 



G(x, t;x,t 



lim[G(x, t;x + r], t, [(f)]) + G{x, t;x-ri, t, [0])] 
2 v^o 



(12) 



we can obtain the following effective action (keeping the same zero energy level as fermions, 
we replace + |A| by i7). 



^ni-q, -n)pR{q, Q) - M'q, -^)PL{q, n) - Vpn{-q, -^)pdq. ^ (13) 
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where T> = Vt^ — — \\\'^ . With the functional expressions (J^) of the fermion fields and 
the effective action (pTSD, we can calculate any correlation functions of the system by taking 
functional average, for example, the fermion Green's function is, 

t- x\ t') = ^j n D<PiDp,G{x, t- x\ t', [0])e-^^['^'''] (14) 

i=R,L 

where Z = J l\i=RLD(piDp,e-'^^'^'Pl As A = 0, we can obtain the usual fermion Green's func- 
tion [iHi. AtV = 0, we can obtain the expectation, < e*'^(*«(^'*)-*i(^''*)) >= (_ --^|L==)"', 
and for a finite V, < e^«(*«(^'*)-*i(^'*)) >~ ( ^ = 1 - which are 



consistent with the exact results [0 where ip = J2/g{^R — 



After integrating out the auxiliary boson fields (f)ji(^L){x,t), we can obtain the collective 
excitation spectrums of the fermion density (A 7^ 0) from the effective action. 



^i = J(l-^(^W + |AP + F(g,A) 



^^2 = J(l-^(^W + |AP-F(g,A) (15) 



where F{q,\) = v (^)^Q'^ ~4|Ap. As V = 0, the energy spectrums reduces to fii 



^2 = \l(f + |Ap; as A = 0, we only have one-branch energy spectrum ^2 = y 1 ~ (2^)^Iq'I) 
which is well-knowen in the usual bosonization representation. Generally, for a finite A, the 
function F(g, A) becomes imaginary for small q (|^| < 1), therefore, the small q excitation 
modes are forbiden. The physics meaning of this phonomenon is unclear, it maybe derives 
from the finite-size and weak-localization of the (anti-) solitons and breathers. 

As A = i-u and including the static tension potential energy of the lattices, the dimer- 
ization parameter u can be determined by taking the minimum of the ground state energy 
of the system [O. When A = Aq and taking Aq as a free parameter, we can obtain its 



renormalization group equation by using the effective action ([13]), and study its low energy 
flow behavior in strong repulsive and attractive interaction regions ||^. On the other hand, 
the renormalized parameter Aqr is related to the mass of the (anti-)solitons and breathers in 
the exact solution of the QSGM. It is clearly shown |^ that the fermion interaction depen- 
dence of the renormalized parameter Ao_r is completely consistent with the exact solution 



of the QSGM. Therefore, we demonstrate that the eigen-functional bosonization method 
is powerful, and the low energy behavior of the strongly correlated systems is controlled 
by the Eikonal-type equations. Even if at one-loop approximation it can give the correct 
results in treating the one-dimensional strongly correlated systems, where in general usual 
perturbation theory is invalid. 

III. CHARGE COLLECTIVE EXCITATION GAP OF THE HUBBARD MODEL 

AT HALF-FILLING 



For the Hubbard model at half filling, it is well-known |[T8| that at half-filling, the spin 
collective excitation gap is zero, and the charge collective excitation gap is not zero, and 
increases with the on-site repulsive Coulomb interaction potential U of electrons. In this 
section, we use the eigen-functional bosonization method to calculate the effective action 
of the electron density fields, so that we can determine the electron interaction dependence 
of the charge collective excitation gap of the Hubbard model at half-filling. This system 
represented by the Hubbard model with electron half-filling has a metal-Mott insulator 
transition at f/ = 0. As U > it is a Mott insulator, and at ?7 = it is a metal. 

In low energy region, and with the linearization of electron spectrum near the Fermi 
levels ^kp, the Hamiltonian of the Hubbard model at half-filling can be written as in the 
continuum limit, 

H = -ikvpYl / dx tlj]i^{x)d^iljRa{x) - iljl^{x)dxtlJLa{x) 

+ V dx [pR^ix) + pL^{x)] [pRlix) + pLlix)] (16) 

where ipRa{x) {ipiaix)) is the right (left) moving electron field with spin a, Pr{l)u{.x) = 
'^R(L)ai^)'^R{L)o-{x) is the electron density field, V = all, and A = all, where a is the lattice 
constant. Here we have used A to present the coefficient of the Umklapp scattering which 
is relevant for repulsive electron Coulomb interaction {U > 0), therefore, in the low energy 
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region renormalized A is different from all, and neglected the backward scattering term 
which is irrelevant for the repulsive interaction V > 0. 

With auxiliary fields (j)Rcr{x,t) and 4>La{x,t) which introduce the constraints pRa{x,t) — 
'^^j^{x,t)'^Rcr{x,t) and pLa{x,t) — '^\^{x,t)'^La{x,t), respectively, and the Hubbard- 
Stratonovich field A_o-(a;,t) = V^Jj^(a;, t)^£,o-(a;, t), the action of the system can be written 



as, 



S = Y, dtd^ {"^lix, t)Ma{x, t)'i/a{x, t) 

c ■' 

+ J dtdxi^j [A^{x,t)A^{x,t) + A*^{x,t)Al{x,t) 
- V \pr^{x, t) + pl^{x, t)] [pRi{x, t) + pLi{x, t)]} 



(17) 



where M^{x,t) = Mo + 0^(x,i), = (V'L(2^> 0> V'L(^> and 

(t)^{x,t) 



A*, Vl 



^ (l)Ra(x,t), A^(x,t)-A ^ 



A*^{x,t)-A\ (j)La{x,t) 



where 'Pi?(L) = idt^idx (choosing h — vp — 1), and A is a parameter (see below). Integrating 
out the electron field ^'^-(a^, ^), we obtain an effective potential, W[(f)] = iTr\n{Ma-) which 
can be represented by using eigen-functionals, 

= iY.t dtdx tr [0.G.(x, t; x\ t', [^0])] | 
G^{x, t; x', t', [0]) = E E T;^y:Li^^ t', m (18) 



where the eigen-functional '^^^i^^{x,t, [(f)]) satisfies the eigen-equation. 



M,(x,t)*il(x,t, 



(19) 



Using Hellmann-Feynman theorem, the eigen- value E^i}^[(j)] reads = ^^^^-^ [0] + E^'^'' [0] , 

where Ei^[0] = dU dtdx^^:il{x,t,[mM^,t)¥:l{x,t,[m, 4?[0] ^ u; T Ek, and 
Ek^^k^ + \A\^. 
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At (pR(L)a-{x,t) = (i.e., V = 0),'we have the solutions corresponding to the eigen- values 
^fc?[0]' respectively, 



{tl 



ikx—iwt 



where a 



(+) 



Uk, -h 



(+) 



.(-) 



Vk,^ = A(^) +zA(^), and tan(2^) = A^^V^^^^- The 



Hubbard-Stratonovich field becomes a constant under this condition, Ao-(x, t) = A, and the 
parameter A is determined by the self-consistent equation, 

A* 



2L^ 



(20) 



k ^A;2 + |A|2 

which only has imaginary solution, A'^'^^ = 0, and A*^'^^ ~ £)g-2nhvp/x^ where D is the band- 
width of the system. The parameter A'^'^^ is a charge gap parameter, and at A = (i.e., 
f/ = 0) it is zero. Therefore, it can be taken as an order parameter of the Mott insulator. 
Generally, we can choose the ansatz. 



(±) 

akuj 



'x,t, [(f)]) = Ak 



(^) 



1/2 



jkx-i{uj+j:l^^'>[4,])t 



(21) 



where ~ 1 is a normalization constant, and the Hubbard-Stratonovich field reads 
A„{x,t) = zA(^)e'^^-'^^^'*^"'^«-'^^^'*\ Substituting (|T]) into the eigen-equation (H), we can 
obtain, 



A«al^) 



q(-) _q(-) q(±)-q(±) . 







(22) 



where Q^^i^iyix, t) = Q%(^q„{x, t) + fl^l^^{x, t). These differential equations can be simplified 
to the Eikonal-type equations. 



(±) 

_Rcr 



-Vr 



(±) 



,A(^)al,^) 



(23) 
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-(A(^))-(e 



(±) 



Z(*.=^e«i*'(P„-P«/<t')(e«^-.'-l) 



where Q° = Q?j, - QL, and Q(±) = Q^^ - Q 

In order to solve these equations, we only keep the linear terms of fji(L)ai^^ ^) and Q\f^^^- 
Under this approximation, we obtain the effective actions Sl^^ [0, p] and 5*1^^ [0, p] of the 
spin and charge parts, respectively. 



s: 



1 



oTT JJ, ■ ' 



c(.) 



(24) 



2 



T j|pijc(.)(g,^^) + PLc(.)(g,^^)h 

where Dels) = l^^-g^-(A(^)±A(^))^ and Xc(5) = XT±Xi where x = (pR{L), Pr{l)- Integrating 
out the auxihary fields 4>rc{s){(1, ^) and (t>Lcis){(l, we obtain the effective actions. 



q,n 

I 



Xe(,)(g,^]) = A, 



c{s) 



V 



V 



V 



4A, 



c{s) 



-q{n + q)T 



V 



(25) 



where p^. 



Pc(s) = {PRc{s),PLc{s)), Ac(s) 



ttD, 



— {q'^ — q^Vl^) ^ With these effective action, we 



(26) 



can obtain the gaps in the spin and charge collective excitation spectrum, respectively, 

A, = 
A, = 2A(^) 

At zero-order approximation (i.e., \^ = 0), the parameter A'-^-' is determined by equation 
(pop, and it increases with the electron interaction strength A (= aU). However, the param- 
eter can be more rigorously determined by taking the minimum of the ground state 
energy of the system, where the ground state energy reads. 
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2(AW)- 
A 



-iTr In (Mq) +-Tr In 



(27) 



where Dc{V) = — (1 + y/7r)g^ — 4(A*^'^))^. By simple calculation, we can obtain, 



IT 



/ 



,/l + V/n 



:ln 



V 



1 + V/nD + ^(1 + y/7r)D2 + (AW)2 
AW 



(28) 



At = 0, it reduces into the equation (PP]). We now can approximately take A = in the 
low energy region, and use equations ( PBD and ( PB| ) to determine the order parameter A*^^) 
of the Mott insulator and the charge collective excitation gap, respectively, where the only 
parameter is the band-width D of electrons. This result is qualititively consistent with the 
exact solution of the Hubbard model at half-filling, where the charge gap increases with the 



on-site Coulomb interaction strength of electrons [18 



IV. PHASE TRANSITION FROM BAND INSULATOR TO MOTT-TYPE 

INSULATOR 

In this section, we study the influence of strong electron correlation on the spin and charge 
collective excitation gaps with a prototype one-dimensional model ||^. This model can be 
used to qualititively explain the origin of ferroelecticity of some transition metal oxides, 
because of a spontaneous lattice dimerization in the strong electron correlation regime. 
By analytical calculations of the spin and charge gaps, we can clearly demonstrate that 
the system has a quantum critical point defined by Ac = at the interaction strength 
V = Vt(Ao), where the order parameter A'^'^^ of Mott insulator and dimerization parameter 
u turns on, and increase with as > Vt(Ao), and the order parameter A*^^-* of band 
insulator goes to a constant as ^ > Vt{Aq). This quantum critical point represents the 
phase transition from band insulator to Mott-type insulator, and as Aq = 0, it reduces into 
the usual metal-Mott insulator transition of the Hubbard model at half-filling. Another 
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quantum critical point defined by = is reached as A'^'^) oo asymptotically for a finite 



Aq. Qualititively, the phenomenal results of Ref. [IS] are consistent with our calculations 



The Hamiltonian of this prototype one-dimensional model reads, 

H = J2{-i (cIq+i. + 4+i,Q.) + Ao(-l)*n,.} + Uj2n,^n,^ (29) 

ia i 

where Cia{cl^) is the electron annihilation (creation) operator with spin cr (=|, |) at site i, and 
njo- = cJ^Cjo- is the electron density operator. The odd and even sites represent oxygen atoms 
(O) and a generic cation (C), respectively, with the energy difference, Ec — Eq = 2Ao. 
It can be easily seen that the Hamiltonian (^) describes different low energy physical 
behavior for the on-site Coulomb interaction U ^ and U oo, respectively. At f/ = 0, 
it represents a simple band insulator where there is an energy gap 2Ao in the electron 
excitation spectrum. However, as f/ ^ Aq, it represents a Mott insulator where charge 
density excitation spectrum has a gap and spin density excitation spectrum is gapless. 
Therefore, there is a transition (mixed- valence) region between the band insulator and Mott- 
type insulator as f/ ~ 2Ao. The finite-size simulations [|,^,^ suggest that the mixed- 
valence region is accompanied by a spontaneous dimerized phase which turns on electron- 
phonon interaction. 

Hep = -u ^(-1)* (cLci+i^ + (30) 



where u is the dimerization parameter. 

In low energy region, and with the linearization of electron spectrum near the Fermi 
levels ifci?, the Hamiltonian (EOl) can be written as in the continuum limit. 



+ V dx [pR^ix) + pL^ix)] [pRiix) + pLiix)] (31) 

-..E/-k.L(.).U.)..L(.)..(.)l 
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where ipRu{^) ii^Laix)) is the right (left) moving electron field with spin a, PR{L)a{x) = 
'^R{L)ai^)'^R(L)a{x) IS the electron density field, V = all, and A = all, where a is the lattice 
constant. Here we have used A to present the coefficient of the Umklapp scattering which 
is relevant for repulsive electron Coulomb interaction {U > 0), therefore, in the low energy 
region renormalized A is different from all, and neglected the backward scattering term 
which is irrelevant for the repulsive interaction V > 0. The spontaneous electron-phonon 
coupling term (|30|) reads in the continuum limit. 



Hep = iuJ2 dx U\^{x)i)La{x) - i)[^{x)i)R^{ 



X 



(32) 



The dimerization parameter u can be determined by taking the minimum of ground state 
energy after including static tension energy of the lattices. 

With auxihary ffelds 4>Ra{x,t) and 0lo-(x, t) which introduce the constraints pRa{x,t) = 
iljRf^{x,t)ipR„{x,t) and pLa{x,t) = 4'Lai^y'l^)i'La{x,t), respectively, and the Hubbard- 
Stratonovich field A_^{x, t) = ip\^^{x, t)'4'L„{x, t), the action of the system without the lattice 
dimerization can be written as, 



S = J2 [ dtdx{^l{x,t)Maix,t)^^{x,t) 

- (pRa{x,t)pRa{x,t) - (f)La{x,t)pLa{x,t)} 
1 



+ J dtdx |- [A|(x,t)A|(x,t) + A*^{x,t)Al{x,t) 

- V [Pi?T ix,t)+ PL] (X, t)] [pRi (X, t) + pLi (X, t)] } 



(33) 



where M^{x,t) = Mq + 0^(x,t), '^U^,t) = (?/;|j^(x, t), ^'[^(x, t)) , and 



Mn 



^ Vr, Ao + a ^ 



Ao + A*, Vl 



^ (j)R.a{x,t), A,{x,t) - A ^ 
^A;(x,t)-A*, 0L,(x,t)^ 



where T^r(l) = idt^id^ (choosing h = vp = 1), and A is a parameter (see below). Integrating 
out the electron field 'if„{x,t), we obtain an effective potential, W[(f)] = iTrln(Mo-) which 
can be represented. 
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^[0] =^T. f'dU dtdx tr {(P^G^ix, t; x', t', [^0])) | 

G.(x, t; x', t', [0]) = E E T;(T7:7T*^^-(^' t, t', [0]) (34) 

where the eigen-functional t, [0]) satisfies the eigen-equation, 

t, [0]) = E^lm^lix, t, [0]) (35) 

Using Hellmann-Feynman theorem, the eigen- value -^^^^^[0] reads -E'f.fc|^[0] = -E'^^'' [0] + S^^'^-' [0] , 
where ^^^"^^[0] = dU dtdx'$';!^il{x,t, [^0])0^(x, t)^2^(x, t, [^0]), independent of u. 
In order to solve the eigen-equation (pSf), we can choose the ansatz, 

^(±)g-i9gQi±'{x,t) 



/I \ 



Vtl/ 



where tan(2^) = A(^V(^o + A^^)), the eigen- value is ^iti[0] = to t Ek + ^^^^[0], E, 
^fc2 + |Ao + AP, = = y'Kl + fc/Efc), -4+^ = = ^i(l-fc/Efc), A, ~ 1 is 
a normalization constant, and A<^(x, t) = Ae'^^--^^'*^~'^«--^'''*\ where A = A^^^ + iA^^^ is 
determined by the self-consistent equation, 

A-'--^E^ (37) 

For a finite Aq, the real part A*^^^ has a non-zero solution (A(-^) < 0), and A^-^) = for 
V < Vt'(Ao), and at Aq = it only has a zero solution A^^^ = 0, where the system reduces 
to the usual Hubbard model at half-filling. Therefore, the parameter A'-^-' characterizes the 
band insulator, and A*^^) characterizes the Mott insulator. The parameters A*^^-* and A*-^) 
can be taken as the order parameters of band insulator and Mott insulator, respectively. 
Substituting (|36D into the eigen-equation (|35|), we can obtain. 



k 



^ A±) ^(^k^^ f o(±'-o(±' .\ A4^^e-*2e / Qi-) _oi-) \ oW-o(±' r. /ooN 
T^rJrJ H tI)- ( e^""^ ~ ) + (±) [e^^- - IJ e'^i- = (38) 



"fc ^ ' o^k 



(-^k f^k 
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where m = [(Aq + A(«))2 + (A«)2]V^ and gg[),(x,t) = Ql^L)A^,t) + fi^l^ix,t). These 
differential equations can be simphfied to the Eikonal-type equations, 



^L^ij — (±r^ — ^(±r /fl<T - — (±) — ^iv,7 + ^o- 

V '^k Pk I \ (^k I 

' ,3(±) (±) \ / (i) \ 

J. J. !m_^ ^o^k ] A±) _ / mo^k ^ A±) ] 75 + 7(±) 



V Pk J ^k 

where = Ql, - Q%, = A/Sj^'^e'^^'e-'^'+f'LV {cQlI-QrI - l)/a^^\ and = 

A*a[?^^e*2^e'3°+4'!^^(e'^B~--"'^^-- - 1)//?^'^^. At A = 0, it is reduced to the same Eikonal-type 
equations as that in the equation (|^). 

In order to solve these Eikonal-type equations, we only keep the linear terms of fj^^yix, t) 
and Q\i^^„- Under this approximation, we obtain the effective actions S'gjj [0, p] and 
S'ljj [0,p] of the spin and charge parts, respectively. 



- \ (pRcis){-q, -^)prc{s) {q, fi) + 4>Lcis){-q, -^)plc{s) {q, ^) 
V 

T -^\pRcis){q,^) + pLc{s)iq,^)\ 



where = - - (m ± A)^, A = ^(A(^))2 + (A«)2 cos {2(6 - 9)) , tan(20) = 

//\(^)^ and Xc(s) = ^Xi where x = <Pr{l), Pr{l)- Integrating out the auxihary fields 
(pRc{s){q,^) and (f)Lc{s){q,^), "we obtain the effective actions, 

1 



Seff.[p] = T^Y. Pl{s){-q^ -^)Xcis){q, ^)pc(s){q, ^) 



q,n 

Xc{s)iq, ^) = ^c(s) 



(41) 



where pjs) = {PRc{s), pLcis)), ^c(s) = -^^(g^ - q^^^) ^ 



With the effective actions 5*1^^ [p] and S^jj-_[p] (^T]), we can obtain the gaps in the spin 
and charge collective excitation spectrums, respectively, 

' ^ ' (42) 

Al = (m + A)' 

At Ao = 0, it reduces to the well-known result of the Hubbard model at half-filling that 
As = and A^ = 4(A'-^^)^, where the charge gap is increased with the electron interaction 
strength U. This can be seen from the self-consistent equation of A^^^ (|37D where A^^^ 
increases with A. Generally, the parameter A can be determined by the minimum of the 
ground state energy, and it is a function of V, X{V). However, we can approximately 
take X = V which does not qualitatively change the behavior of the spin and charge gaps 
determined by equations (|37| ) and (P2|), and further simplifies our calculation. Equation ( ^21) 
clearly shows that the system has two quantum critical points defined by the charge and 
spin gaps, one is at A^ = 0, and another is at A^ = 0. The physical behavior around the 
quantum critical points can be understood by solving the self-consistent equation (^). By 
simple calculation, we can obtain the following solutions, 

0, 



^inhvp/V — 1 I 4 

[a(^)(\/), V<Vt{A, 



A 



1/2 



V < Vt{Ao) 

V > Vt{Ao) 



1 



V > Vr(Ao) 



(43) 



where Vt(Ao) = 2TchvF/\n[2{D + ^jD'^ + Ag/4)/Ao], and A^^'>{V) is determined by equation 
(^) with A^^) = 0. The interaction V = Vt{Aq) corresponds to the quantum critical point 
Ac = 0. It means that there exists a phase transition at Ac = where the parameter Ac 
shows different behavior on the two sides of the point V = Vt{Aq). As < Vt{Aq), the 
parameter is zero, and A'^'^^ decreases with increasing V. In this case, the parameter 
A is real. On the other hand, as ^ > Vt(Ao), the parameter A'^'^^ becomes a constant, and 
A*^'^'' is not zero and increases with V. The parameter A becomes a complex quantity. For 
the strong interaction V > Vt(Ao), the charge and spin gaps can be written as, respectively. 
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^ 4(A">)^ 



^0 



It is noted that for weack interaction V < Vj'(Ao) the spin gap retain invariant, and the 
quantum critical point = corresponds to A'^'^^ oo (i.e., V, D oo). In fact, this 



is due to our simple approximations that the self-consistent equation (|37|) is obtained at 
0o-(x,t) = (i.e., V = 0), it does not include the high-order correction of V. However, our 
results are qualitatively correct after including this influence. 

Including the electron-phonon interaction, the equation (^) keeps invariant, but the 
parameters m, 9 and the self-consistent equation of A are changed as, 

m= [(Ao + A(^))2 + (A«-m)Y/', tan(2^^) 



Ao + A(«) 



A ^ Ao + A(^) A ^ A(^) - u 



where Ek = [k^ + m2]V^ A = |A| cos (2{9 - ^)), and tan(2^) = A^^V^^'^^- Here we only 
consider the static uniform dimerization of the lattices. Including the static tension energy 
of the lattices ^au"^, the dimerization parameter u can be determined by the self-consistent 
equation, 

2 _ A(^) + A /A(^) -M \ ^ / 1 1 \ 

where z = |A| sin (2{e - 6)) cos2(2^)/(Ao + A^^)), E<'\V) = [(1 ± V/7i)q^ + A^^^^J^s, and 
E^^'^^ = [g^ + A^^^j]^/^. This equation is derived from the minimum of the ground state 
energy, duEg[u] = 0. According to equations ( ^5] ) and (^Uf), we see that the dimerization 
parameter u, similar to A^^^ is zero for weack interaction V < Vt{Aq), and turns on at 
V = Vt(Ao). However, for strong interaction V > Vt(Ao), the parameter A*^^^ basically 
keeps a constant, because the dimerization parameter m is a small quantity, and it goes to 
zero as V > Vt(Ao). 
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In usual transition metal oxides, there is the strong hybridization between d-orbit elec- 
trons of cation and p-orbit electrons of oxygen atom. Therefore, at enough strong elec- 
tron correlation, there spontaneously takes place the lattice dimerization and deformation- 
induced charge transfer, which is possibly the origin of ferroelectricity of the transition metal 
oxides. Our calculations demonstrate that in the strong interaction region V > Vt{Aq), the 
system has the spontaneous lattice dimerization, which is consistent with the previous nu- 
merical calculations |]8|,p0|,^ . 



V. DISCUSSION AND CONCLUSION 

In general, for a one-dimensional strongly correlated electron system, if there is only 
low energy excitation modes with small momentum near the Fermi levels ±kp, it shows 
the Luttinger liquid behavior. If there are disorder impurities which induces the new low 
energy excitation modes with large momentum transfer between two Fermi levels ±kp, the 
electrons are localized, and the system becomes an insulator. At electron half-filling, there 
appears the Umklapp scattering which represents the low energy excitation modes with large 
momentum 4k p transfer between the Fermi levels ±kp, the system becomes a Mott insulator. 
Therefore, the low energy excitation with large momentum transfer between the Fermi levels 
±kp (i.e., crossover excitation) makes the system become the insulator for repulsive electron 
interaction. 

Usual perturbation methods fail to treat such the system, because the crossover excitation 
term(s) is relevant, and determines the low energy behavior of the system. It is very desirable 
to find a new method which can exactly and effectively treat this kind of systems. The 
eigen-functional bosonization is a good candidate, because it has more advantages than 
usual perturbation theory: a). Without the crossover excitation term(s), just as usual 
bosonization methods it can exactly treat the system, b). For an interacting-free 

electron system, it is exact for dealing with the crossover excitation term(s). c). With this 
method, the problem of the strongly correlated electron systems ends in to solve the Eikonal- 
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type equations, which can be exactly solved by a series expansion of the boson fields (5° (x, t), 
and/or by a computer calculations, d). For some integrable systems, such as the quantum 
sine-Gordon model and the Hubbard model at half-filling, at one-loop approximation (i.e., 
only keeping the linear terms of the Eikonal-type equations) , it can give the correct results 
that are qualititively consistent with the exact solution of the systems for both the weak 
and strong electron interactions. 

It is the key point of the eigen-functional bosonization method to solve the Eikonal-type 
equations of the eigen-functional, and use these eigen-functionals to calculate the electron 
Green's functional. With the electron Green's functional, we can calculate the ground state 
energy and the effective action of the system as well as correlation functions. This method 
can also be used to treat one-dimensional electron-phonon interaction systems, two cou- 
pled spin-chain or quantum-wire systems, and one-dimensional disorder interacting electron 
systems. It can also be extended to two- and three-dimensional electron systems. 
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